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Abstract 

We study the deterministic and stochastic rotational dynamics of ferromagnetic nanoparticles in a 
precessing magnetic field. Our approach is based on the system of effective Langevin equations and 
on the corresponding Fokker-Planck equation. Two key characteristics of the rotational dynamics, 
the average angular frequency of precession of nanoparticles and their average magnetization, are 
of our interest. Using the Langevin and Fokker-Planck equations, we calculate both analytically 
and numerically these characteristics in the deterministic and stochastic cases, determine their 
dependence on the model parameters, and analyze in detail the role of thermal fluctuations. 
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I. INTRODUCTION 


Ferromagnetic single-domain nanoparticles possess a number of unique properties, such, 
for example, as superparamagnetism [1-3], giant magnetoresistance [4, 5] and quantum tun¬ 
nelling of magnetization [6-8]. These and other properties provide a basis for numerous 
current and potential applications of magnetic nanoparticles in data storage [9-11], spin- 
tronics [12, 13], drug delivery [14-16] and hyperthermia [17-20], to name only a few. The 
magnetization dynamics plays a central role in most of these applications, and its charac¬ 
teristics strongly depend on whether the nanoparticles move or not. In the latter case, the 
time evolution of magnetization can often be described phenomenologically by the Landau- 
Lifshitz or Landau-Lifshitz-Gilbert equation [21, 22], By adding the thermal torque [3] and 
spin-transfer torque [23, 24] to these equations, they can also be used to study thermal and 
spin-transfer effects. Within this framework, a wide variety of phenomena, including preces- 
sional switching, self-oscillations and thermal relaxation of the nanoparticle magnetization, 
have already been investigated (for a review, see Ref. [25] and references therein). 

Although there is experimental evidence that ferromagnetic nanoparticles can freely ro¬ 
tate even in a solid matrix [26], it is more obvious that the former case occurs for nanopar¬ 
ticles suspended in a fluid. These systems, also known as ferrofluids, exhibit a number of 
unique properties [27, 28]. In dilute suspensions, some of these properties are completely 
determined by the magnetic and mechanical dynamics of independent nanoparticles. Re¬ 
markably, in the case of high-anisotropy nanoparticles the magnetization motion describes 
the nanoparticle rotation as well. Because of its simplicity and efficiency, this approach is 
very useful for investigating thermal effects in these systems (see, e.g., Refs. [29, 30]). In 
particular, it has been used to predict and study the thermal ratchet effects in ferroflu¬ 
ids subjected to a linearly polarized magnetic field [31, 32], and to determine the specific 
absorption rates [33, 34], 

In this work, we present a detailed study of the rotational dynamics of highly anisotropic 
nanoparticles in a precessing magnetic field. We focus on the average angular frequency of 
nanoparticle rotation and on the average nanoparticle magnetization. The dependencies of 
these quantities on the model parameters, especially those that exhibit qualitatively different 
behavior with and without thermal fluctuations, are our main interest here. 

The paper is organized as follows. In Sec. II, we describe the model and main approxima- 
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tions and write the basic system of Langevin equations governing the rotational dynamics 
of ferromagnetic nanoparticles. The corresponding Fokker-Planck equation and the system 
of effective Langevin equations, which is more convenient than the basic one, are obtained 
in Sec. III. In the same section, we show that the interpretation of multiplicative Gaussian 
white noises in effective Langevin equations does not influence the statistical properties of 
their solution. Section IV is devoted to the computation of the average angular frequency of 
precession of nanoparticles and their average magnetization in the deterministic case. The 
effects of thermal fluctuations are considered in Sec. V. Here, we confirm the validity of the 
system of effective Langevin equations, use it for calculating the above mentioned average 
characteristics, and discuss the role of thermal fluctuations. Finally, our main conclusions 
are summarized in Sec. VI. 

II. MODEL AND BASIC EQUATIONS 

We consider a spherical ferromagnetic particle of radius a, which rotates in a viscous fluid 
under a uniform magnetic held H = H(f). In our study we use the following assumptions. 
First, the exchange interaction between magnetic atoms is assumed to be so large that 
the magnitude |M| = M of the particle magnetization M can be considered as a constant 
parameter. Second, the particle radius is assumed to be so small (less than a few tens 
of nanometers) that the nonuniform distribution of magnetization becomes energetically 
unfavorable, i.e., a single-domain state with M = M(f) is realized. And third, the magnetic 
anisotropy held is assumed to be so strong that the magnetization is directed along this held, 
implying that M is frozen into the particle body. With these assumptions, the rotational 
dynamics of a ferromagnetic particle is governed by a pair of coupled equations 


M = cu x M, 

(2.1a) 

Jib = VM x H — 6//Utu. 

(2.1b) 


Here, cu = cu(f) is the angular velocity of the particle, the overdot denotes the time deriva¬ 
tive, J = (2/5 )pVa 2 is the moment of inertia of the particle, p is the particle density, 
V = (4/3)7ra 3 is the particle volume (we associate the hydrodynamic volume of the particle 
with its own volume), r/ is the dynamic viscosity of the fluid, and the cross denotes the 
vector product. 
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The former equation in (2.1) is a special case of the kinematic relation a = cu x a, which 
holds for an arbitrary frozen vector a of a fixed length, and the latter one is Newton’s second 
law for rotational motion. The first and second terms in the right-hand side of Eq. (2.1b) 
are the torques generated by the external magnetic field and viscous fluid at small Reynolds 
number (< 2 x 10 3 ), respectively. Because the particle size is sufficiently small, the left-hand 
side of this equation, i.e., the rate of angular momentum J cu, can safely be neglected in a 
wide frequency domain. Using this massless approximation and assuming that a random 
torque £, = £,(£), which is generated by the thermal motion of fluid molecules, is also applied 
to a nanoparticle, we obtain 


cu = —M x H -|- 

br) 6r/V 


( 2 , 2 ) 


With this result, Eq. (2.1a) reduces to the equation (see, e.g., Ref. [30]) 


M = ~ M x (M x H) - -L M x £,, 
br) br/V 


(2,3) 


which describes the stochastic rotation of ferromagnetic nanoparticles in a viscous fluid. Note 
that, in spite of the similarity in appearance, Eq. (2.3) strongly differs from the stochastic 
Landau-Lifshitz equation describing the magnetization dynamics of fixed nanoparticles. The 
main difference leading to a qualitatively different behavior of M is that Eq. (2.3) does not 
contain the gyro magnetic term in the deterministic limit. In particular, it is this term that 
is responsible for the magnetization of nanoparticle systems in a rotating magnetic field [35]. 

Since the nanoparticle magnetization M does not depend on time, it is convenient to 
rewrite Eq. (2.3) in spherical coordinates. To this end, we first represent the magnetization 
vector as M = Mm with 


m = sin 9{ cos (p e x + sin e y ) + cos 9 e z , (2.4) 

where 9 = 9{t) and p = ip(t) are the polar and azimuthal angles of the nanoparticle mag¬ 
netization, respectively, and e x , e y and e, are the unit vectors along the corresponding axes 
of the Cartesian coordinate system xyz , whose origin is located at the nanoparticle center. 
Then, introducing the rescaled random torque as 

/ 1 \ V 2 

C '~ \12r/Vk B Tj ^ 


4 



where ks is the Boltzmann constant and T is the absolute temperature, from Eq. (2.3) one 
can obtain the following basic system of stochastic Langevin equations: 


ft ±dw 12 

0 = -— - \j — (Cx Sin p - Cy COS p) , 


Ti 89 v t 2 

1 dw 




^“[(Cx COS p + Cy sin p) COt 9 Cz] ■ 


( 2 . 6 ) 


t i sin 2 9 dp 

Here, w = W/M 2 , W = -M • H is the Zeeman energy density, the dot denotes the scalar 
product, and 


T\ = 


6r/ 


t~2 = 


6r]V 


(2.7) 


M 2 ' ‘ k„T 

are the time scales characterizing the nanoparticle rotation induced by the external magnetic 
held and thermal torque, respectively (t 2 /2 is also called the Brownian relaxation time). The 
Cartesian components ( u {v = x, y, z ) of C are assumed to be independent Gaussian white 
noises with zero mean, (£„) = 0, and correlation function (Cv(t)Cv(t')) = A S(t — t'), where 
(•) denotes averaging over all realizations of Wiener processes W„(t) producing noises ( u (for 
more details, see the next section), A is the dimensionless noise intensity, and 5(t) is the 
Dirac 5 function. 

Finally, we choose the external magnetic held in the form 


H = H[cos(cot)e x + sin(a;i)e 3/ ] + H z e z , (2.8) 

where H and u are the amplitude and angular frequency of the circularly polarized (rotating) 
component of H. and H z e z is the constant component of H (see Fig. 1). In this so-called 
precessing magnetic held, the reduced energy density w = — m • h (h = H jM) is written as 


w = —h sin 9 cos(uit — p) — h z cos 9 
with h = H/M > 0 and h z = H z /M > 0. 


(2.9) 


III. FOKKER-PLANCK EQUATION 

An important feature of the Langevin equations (2.6) is that the noises (u are multi¬ 
plicative, i.e., they are multiplied by functions of the angles 9 and p. It is well known (see, 
e.g., Refs. [36, 37]) that the statistical properties of one-dimensional systems described by 
Langevin equations with multiplicative noises depend on the noises interpretation. In con¬ 
trast, the statistical properties of some multi-dimensional systems do not depend on how the 
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FIG. 1. (Color online) Schematic representation of the model. Nanoparticle of radius a with frozen 
magnetization M rotates about the origin of the Cartesian coordinate system in a viscous fluid 
under the action of a precessing magnetic field H. 


multiplicative noises are interpreted [38]. Therefore, to determine if the noises interpretation 
influences the statistical properties of 9 and ip and to And the Fokker-Planck equation for 
the probability density of these angles, Eqs. (2.6) must be specified more precisely. 

For this purpose it is convenient to rewrite the system of stochastic equations (2.6) in the 
form 

3 

Ui = fi(u,t) + 5Zsij(u)Cj. (3- 1 ) 

3 = 1 


Here, iq {i = 1, 2) are the elements of the 2x1 matrix [two-component column vector 
u = u(f)] ( Ui ) = (ul) with u\ = 9 and u 2 = <p, the drift terms fi(u, t) are the elements of 
the 2x1 matrix 


Ui) 


* dw/dui ^ 
^sin -2 iq dw/du 2 j 


(3.2) 


with w taken from equation (2.9) in which the angles 9 and (p are replaced by the variables 
Mi and u 2 , respectively, Ci = C xif), C2 = CC3 — C z(t), and the functions ^(u) are the 
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elements of the 2x3 matrix 


( 9ij ) 



( 

\ 


sin U 2 — cos 112 

cot Mi cos u 2 cot Mi sin u 2 


0 N 

-V 


(3.3) 


Then, to specify Eqs. (3.1), we first assume that the increments 5uj = Ui(t + r ) — Ui(t ) of 
the variables iq at r <C min{ri, T 2 } are given by 

3 

5ui = fi( u, t)r + ^2 9ij[ u (t + A jT)\SWj, (3.4) 

3 =1 


where A j € [0,1] are the parameters characterizing the action of white noises Q, and SWj = 
Wjit + r) —Wj(t) are the increments of Wiener processes Wj(t) generating Q. Because these 
noises are assumed to be independent and statistically equivalent, the increments SWj can 
be completely characterized by two conditions 


(SWj) = 0, {SWjSWi) = A SjiT 


(3.5) 


with Sji being the Kronecker delta. Finally, taking into account that Uk(t + \jt) 
XjSuk and expanding the last term in Eq. (3.4) to linear order in r, we obtain 

3 


Sui 


fi(u,t)r + ^2g ij (u)5W j 
3 =1 




%^^,(u )SW l SW j . 


~ u k (t) + 


(3.6) 


Thus, the stochastic equations (3.1) are specified by the difference scheme (3.6) in which 
the noises action is accounted for not only through the increments SWj of Wiener processes 
generating Q, but also by the parameters A j realizing an addition connection of the system 
with these white noises. Since the last term in the right-hand side of Eq. (3.6) is of the order 
of r [cf. Eq. (3.5)], this connection is able to strongly modify the statistical characteristics of 
Ui. Although the cases with A j = 0,1/2, and 1 that correspond to the Ito [39], Stratonovich 
[40], and Klimontovich[41] interpretations of Langevin equations, respectively, are usually 
considered, any other values of A j are allowed from a mathematical point of view. Therefore, 
the choice of the parameters A j for Eqs. (2.6) can only be made on physical grounds (see 
below). 
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Now, using Eqs. (3.5) and (3.6) and the two-stage procedure of averaging [42], we can 
derive the Fokker-Planck equation that corresponds to the Lange vin equations (3.1). Intro¬ 
ducing the probability density P = P(U, t) that u(f) = U as P = (6 (u(f) — U)), where U 
is a constant column vector with components U\ and t/ 2 , the straightforward calculations 
[38] lead to the following Fokker-Planck equation: 

+ + M uj)p 


i= 1 


A 2 3 

-Iee 


d 2 


2 "h ^ du * du * 

t,k= 1 j =1 


»y(U)9/y(U)-P = 0, 


(3.7) 


where 


M U) = A EE a 


3gy(U) 

0U t 


9kj(U) 


(3,8) 

(IIJu 

k=l j =1 

are the additional noise-induced drift terms that depend on the interpretation (i.e., values 
of the parameters A j) of stochastic equations (3.1). It should be noted, however, that since 
the noise (3 is additive and so dg^{\J)/dUk = 0, these terms and, as a consequence, the 
probability density P do not depend on A 3 . 

If the reduced magnetic energy w does not depend on time, then /j(U, t) = /*(U) and P 
tends to the equilibrium probability density Po = Po(U) as t —> oo. In this limit, Eq. (3.7) 
for P 0 reads 

d 


i= 1 


S^( /l(U) + /l(u) ) Pc 


A A d 2 


?EE 


- gij (XJ)g kj (XJ)P 0 = 0. 


(3.9) 


2 ^ dU,dU k 

i,k =1 j =1 * K 

It is natural to assume that the solution of this equation is the Boltzmann probability 
density, which for h = 0 can be written in the well-known form 


nh z 


Po = T • \ sin U\ e KhzCOsUl , (3.10) 

47T smh(Kh 2 ) 

where k = t 2 /ti = M 2 V/(k B T). Substituting Eq. (3.10) into Eq. (3.9) and using the 
definitions (3.2) and (3.3), we straightforwardly obtain 

kvh 

Ai + A 2 — 1 H—t - [ 2 + A(Ai + A 2 — 3)] cos LA 


A 

(A 2 — Ai) cos(2f/ 2 

( K h ' 2 


sin 2 IJ { 


+ nh z cos Ui 


A 


(A - 1) sin 2 U k = 0. 


(3.11) 



This condition holds for all possible values of the variables U\ and P 2 (0 < U\ < tt, 0 < P 2 < 
2n) and parameter nh z (0 < nh z < 00), i.e., Eq. (3.10) is the solution of the Fokker-Planck 
equation (3.9), only if 

A = 1, Ai = A 2 = —. (3.12) 

Thus, if Eqs. (2.6) with Gaussian white noises of unit intensity are interpreted in the 
Stratonovich sense, the random rotations of nanoparticles are characterized by Boltzmann 
statistics at long times. 

Now, using the conditions (3.12) and introducing the variables © = U\ and <f> = C/ 2 , the 

Fokker-Planck equation (3.7) can be rewritten in the form 

<9P 1 <9 f dw 1 A 1 <9 / dw A 

dt t 1 (90 y <90 k C °^ ~ J t\ sin 2 0 <9<f> \ (9$ / 


(3.13) 


1 ( <9 2 1 «9 2 A „ 

V 2 \d& + sin 2 0 d§ 2 ) P ~ °' (3 ' ] 

We assume that the solution P = P(@, <£>, i) of this equation is properly normalized, i.e., 


<90 / d$P(0,$,f) = 1, 


(3.14) 


and satisfies the initial condition P(@,T,0) = h(0 — ©o)h( < h — $ 0 ) with @ 0 = 9(0) and 

$0 = <p(0). 


A. Effective Langevin equations 

According to the above results, the basic Langevin equations (2.6) should be interpreted 
in the Stratonovich sense. Due to this fact and because the system of two equations (2.6) 
contains three Gaussian white noises, the study of the rotational dynamics of nanoparticles 
by the numerical solution of these equations is not quite practical. Therefore, it is convenient 
to use, instead of Eqs. (2.6), a system of effective Langevin equations satisfying the following 
requirements. First, the statistical properties of solutions of the basic and effective equations 
must be the same and, second, the effective equations must be interpreted in the Ito sense 
and contain two rather than three independent Gaussian white noises, ft has been shown 
[43] that the corresponding system of effective Langevin equations can be written as 

• 1 dw 1 ^ P2 

0 — — H cot 0 + \ — Hi, 

Ti 06 r 2 V T 2 

(3.15) 

1 dw [2 1 

<P ~ -ryyv-b \ -:—^ A*2, 

Ti sm 6 OLp V A sm t) 
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where /i; = (i = 1,2) are independent Gaussian white noises with zero means, 

{(J>i(t)) = 0, and delta correlation functions, Hi(t')) = S(t — t'). Note that the similar 
system of effective Langevin equations, which corresponds to the Landau-Lifshitz-Gilbert 
equation describing the stochastic dynamics of magnetization in single-domain ferromag¬ 
netic nanoparticles embedded into a solid matrix, has been proposed in Ref. [44], 

According to the results of Ref. [43], the probability density of the solution of Eqs. (3.15) 
satisfies the Fokker-Planck equation (3.13). As a consequence, the rotational properties 
of ferromagnetic nanoparticles can be described either by Eqs. (2.6) interpreted in the 
Stratonovich sense or, equivalently, by Eqs. (3.15) interpreted in the Ito sense. A remarkable 
feature of the latter equations is that, independent of their interpretation, the corresponding 
Fokker-Planck equation is given by Eq. (3.13). Indeed, rewriting Eqs. (3.15) in the form 

2 

u t = Fi( u, t) + y G i:j (u)nj, (3.16) 

l=i 


where 


(A) = - 


Tl 


\ 


dw/dui — (1/ft) cot 
^ sin -2 Ui dw / du ,2 J 


and 


(Gy) = 


A 


\ 


0 1/ sin u i 


one can straightforwardly verify that the condition 




k=l 


(3.17) 


(3.18) 


(3.19) 


holds for all i and j. This means [38] [see also Eqs. (3.7) and (3.8)] that the Fokker-Planck 
equation associated with Eqs. (3.15), i.e., Eq. (3.13), does not depend on the parameters A j 
providing a quantitative interpretation of stochastic equations (3.15). While this conclusion 
is obvious for Ai (because the noise is additive), the independence of Eq. (3.13) on 
A 2 is rather surprising (because the noise yU 2 is multiplicative). We note in this context 
that, in contrast to the univariate case, there always exists a class of multivariate Langevin 
equations with multiplicative Gaussian white noises whose interpretation does not influence 
the corresponding Fokker-Planck equations [38]. The above results show that Eqs. (3.15) 
belong to this unique class of Langevin equations. 


10 



IV. NOISELESS CASE 


Before proceeding with the study of thermal effects, we first briefly discuss the deter¬ 
ministic (noiseless) case. In this case, taking the limit T 2 —> 00 and using Eq. (2.9), from 
Eqs. (2.6) we obtain the following system of deterministic equations: 


T[0 = h cos 9 cos t/j — h z sin 9, 
sin tj} 


(4.1) 


Ti^p — Q — h 


sin 9 ’ 


where ^ = ut — (p is the lag angle and fl = uti is the dimensionless frequency of the rotating 
component of the magnetic field (2.8). Assuming that 9 —)■ 9^ = const and if} —> if}^ = const 
as t —> 00 , the above differential equations in the long-time limit are reduced to a set of 
trigonometric equations 


h cos 6*00 cos ^oo — h z sin 9^ = 0 , 

12 sin doo — h sin i/j^ = 0 . 

Because the solutions of these equations in the cases with h z > 0 and h z 
different, we consider them separately. 


(4.2) 


0 can be quite 


h z > 0 


Using Eqs. (4.2) and the condition h z > 0, it can easily be shown that the stationary 
solution of the deterministic equations (4.1) is given by 


(9 \ 


\^Poo J 


= arcsm 


— \/r 2 - Vr 4 - 4f w 

V2 


UvJ 


(4.3) 


(r 2 = fl 2 + h 2 + hi). It can also be proven that this solution is stable with respect to small 
perturbations of angles 9 and Therefore, the solution of Eqs. (4.1) at h z > 0 always tends 
to the stationary solution (4.3) as t —> 00 . In particular, the ^-component of the reduced 
nanoparticle magnetization in the long-time limit, m z = cos 9^, can always be represented 
in the form 

1 / ~ 

(4.4) 


m z = 




20 2 - r 2 + Vr 4 - 4Q 2 /i 2 . 


In general, m z as a function of the parameters U, h and h z exhibits the expected limiting 
behavior: rn z —> h z /\Jh? + h 2 z as Vt —>■ 0, rn z —> 1 as fl —> 00 , h —> 0 or h z —> 00 , and 
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FIG. 2. (Color online) Nanoparticle magnetization /x and nanoparticle angular frequency x as 
functions of the reduced magnetic field frequency £l/h. The solid curves represent the analytical 
results (4.5) and (4.7), and the symbols indicate the numerical results obtained at t = 10 7 ti by 
solving Eqs. (4.1) for maghemite nanoparticles in water. The numerical values of /x and x-> which 
are calculated at h z = 10~ 3 and h z = 0, respectively, do not depend on the initial values 0(0) and 
-0(0) of the polar and lag angles. 


m z —)■ 0 as h —> oo. But the dependence of m z on hi and h under the condition that h z —> 0 
is not so obvious. Indeed, from Eq. (4.4) one obtains 


0, Q/h< 1 

/x = lim m z = 

v /i - a/h > i, 


(4.5) 


i.e., the nanoparticle magnetization /i depends only on the ratio Q/h and, what is more 
important, the behavior of /x in the regions fl/h < 1 and fl/h > 1 is qualitative different. As 
illustrated in Fig. 2, the numerical solution of Eqs. (4.1) obtained by the fourth-order Runge- 
Kutta method confirms this theoretical result. Here and in the following, the numerical 
calculations are carried out for maghemite ( 7 -Fe 2 C> 3 ) nanoparticles in water. In this case, 
the saturation magnetization of nanoparticles, dynamic viscosity of water and characteristic 
time t i at room temperature T = 298 K are given by 47 tM = 3.89 x 10 3 G, = 8.90 x 10~ 3 P 
and Ti = 5.56 x 10~' s, respectively. 
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B. h z = 0 


Equation (4.5) shows that the case with h z 
Eq. (4.3) yield 


0 is special. If Q/h < 1 then Eqs. (4.2) or 


9 


OO 


7r 
2 ’ 


Ipoo 


= arcsin 



(4.6) 


and, as in the previous case, this solution is stable. In contrast, if Vt/h > 1 then the steady- 
state solution of Eqs. (4.1) is periodic in time with period t st = nri/y/fi, 2 — h 2 (periodic 
regime of rotation) [45]. More precisely, in this case the angles 9{t) and ij)(t) are changed 
in such a way that 9(t + t st ) = 9{t) and ij;(t + t st ) = n + if}it). Using these results, it is 
possible to determine the reduced angular frequency of nanoparticles, which is defined as 
X = (1 /oj) lirry^oo c p(t)/t. Indeed, since (p(t) = cut — ij;(t), from this definition one gets x — 1 
for f1 /h < 1 and x = 1 — vr /(cot st ) for Q/h > 1 , i.e., 


1 , Q/h < 1 

X = { ,_ (4.7) 

[i - yi - (KJW, a/h > i. 

This dependence of x on Vt/h is also in excellent agreement with the numerical results, as 
shown in Fig. 2. 

Comparing Eq. (4.7) with Eq. (4.5), we can see that the nanoparticle magnetization // and 
the nanoparticle angular frequency x are connected in a remarkably simple way: /x + x = 1 . 
It should be noted that, although the steady-state dynamics of the unit magnetization 
vector m at h z = 0 and Q/h>l may strongly depend on the initial direction of this vector 
(see Fig. 3 for illustration), there is no initial-state dependence for x- Thus, the condition 
H+X — 1 is universal in the sense that it holds for all possible values of the reduced magnetic 
field frequency VL/h and does not depend on the initial values 9(0) and ^(0) of the polar 
and lag angles. 

To avoid confusion in interpreting the above condition, we first recall that at h z > 0 the 
stationary solution (4.3) of Eqs. (4.1) is stable for all values of the ratio f 1/h. In this case, 
m z is given by Eq. (4.4), x = 1 and, in addition, m z can be approximated by Eq. (4.5) 
if h z is small but nonzero. In contrast, since at h z = 0 the stationary solution (4.6) of 
Eqs. (4.1) is stable only if Q/h < 1 and these equations at Q/h > 1 have a periodic steady- 
state solution, in this case the nanoparticle angular frequency is given by Eq. (4.7) and the 
nanoparticle magnetization equals zero: n — 0. The last result following from the definition 
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FIG. 3. (Color online) Examples of steady-state trajectories of the rotational motion of nanopar¬ 
ticles illustrating their dependence on the initial polar angle 0(0). The trajectories are obtained 
from the numerical solution of Eqs. (4.1) for = 2.5, h = 1, h z = 0, ip(0) = 0, and 0(0) = 0.4 (a), 
0(0) = 7 t/2 ( b ), 0(0) = 2.7 (c). In this figure, the angles are measured in radians, the lag angle 
-*/>( t ) is reduced to the interval [0,7r], the arrows indicate the direction of time evolution, and the 
degenerate trajectory b corresponds to the nanoparticle rotation in the xy-plane [i.e., 0(t) = vr/2]. 

/i = (l/f st ) /g St dt cos 0(f), which accounts for the existence of periodic solution of Eqs. (4.1) 
at fl/h > 1, shows that the rotating magnetic field (when h z = 0) does not magnetize the 
reference systems. Thus, the condition /x + y = 1 holds if /i is associated with m z at h z <C 1 
(not at h z = 0) and y is taken at h z = 0. We note in this context that the rotational regime 
of nanoparticles, which exists at infinitesimally small h z , is completely destroyed by thermal 
fluctuations (see below). 

V. EFFECTS OF THERMAL FLUCTUATIONS 

Next, to study thermal effects in the rotational dynamics of ferromagnetic nanoparticles 
with frozen magnetization, we solve analytically the Fokker-Planck equation (3.13) and 
numerically the system of effective Langevin equations (3.15). 

A. Steady-state solution of the Fokker-Planck equation 

The results obtained for the noiseless case suggest that, depending on the model param¬ 
eters, the steady-state solution P st of the Fokker-Planck equation (3.13) at t —> oo can be 
represented as a function of two variables 0 and T = cut — <h, i.e., P st = P s t(0, T). Using the 
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relation dP st /dt = u>dP s t /<9\l/ and Eq. (3.13), we can find the equation for the steady-state 
probability density P s t directly from Eq. (3.13). For brevity, it is convenient to write this 
equation in the operator form 


LRt = hifl 


dP 


st 


<9 t ’ 


(5.1) 


where the Fokker-Planck operator L is defined as 


f <9/(9 w \ k d f dw 

LPs t = 7WT K7^: - cot 0 P st + . 2 

<90 V 90 / sin 0 d^f \ d^f 


02 1 d 2 |p 

I . O ^ r\.T.o I St 


(5.2) 


<9@ 2 ' sin 2 0 (9'k 2 

with to = w(0,tf) = —h sin 0 cos T — h z cosQ. In particular, if O = 0, then P st is reduced 
to the equilibrium Boltzmann probability density 


P 0 = P o (0, ) = - sin 0 e - ' tTO ( e >’® ,) 

Z 


(5.3) 


(Z = d(~) d if sin 0 e which is the normalized solution of the equation LP 0 = 

0. 

Assuming that kQ <C 1, the steady-state probability density P st can be expanded in a 
power series of kQ. In the linear approximation in nil this expansion yields 


Pt — (1 + kQF)Pq, 

where, according to Eq. (5.1), F = P(@, \k) is the solution of the following equation: 

dP 0 


L(PqF) = 


<9T 


(5.4) 


(5.5) 


Since the probability densities P st and Pq are normalized, the function P must also satisfy 
the condition 

/»7T /*27T 

dO d'hP(0, \H)Po(0, \k) = 0. (5.6) 

Jo Jo 

In what follows, we restrict ourselves to the case when h z = 0 and nh <C 1. Then, using 
Eqs. (5.3) and (5.5), it is not difficult to show that in the main approximation in nh the 
function P is determined by the equation 


1 d ( . ^dF\ 1 8 2 F , . ^ . T 

Sin0— I -|— o „ tt-t-t = — nh sm 0 sin w. 


sin 0 dG V ^ ” <90 J ' sin 2 0 <9T 2 
The solution of this equation, which vanishes as nh —> 0, has the form 

P = —nh sin 0 sin T. 


(5.7) 


(5.8) 
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Therefore, taking into account that, up to quadratic order in K,h, Z = 47t(1 + n 2 h 2 /6) and 


Pn = 


sin 0 


K' 


■h 2 


47t L 

from Eq. (5.4) one immediately gets 


1 + nh sin 0 cos T — (1 — 3 sin 2 0 cos 2 T) 

6 


(5.9) 


P st = P 0 H- n 2 hVt sin 2 0 sin T. 

87T 


(5.10) 


To avoid any confusion, we emphasize that this result is obtained under the assumption that 
P s t = P s t(@, d'), h z = 0 and n max (fl, h ) <C 1. 


B. Simulation results 


Introducing the dimensionless time t — t/r\ and using Eq. (2.9), the system of effective 
Langevin equations (3.15) in the rotating frame can be written as 


d9 

dl 


h cos 9 cos ip — h z sin 6 + —cot 6 + 

hi 



(5.11) 


dip sin ip / 2 1 

If = “ - h ~ 7) ~ \ ~ d2, 

dt sm 9 V n sm cl 

where = /q(t) = y/riHiitTi) (i = 1,2) are dimensionless Gaussian white noises with 
(jdi(t)) = 0 and (/b(f) jdi(t')) = 5(t — t'). It is this system of equations that we used in our 
simulations. 

In order to verify if thermal fluctuations are properly taken into account in the effec¬ 
tive Langevin equations (5.11), we solved these equations by the Runge-Kutta method and 
calculated the quantity 


a = 



hi~n 


(5.12) 


According to Eq. (5.10), this quantity characterizes the difference between the steady-state 
(when Q / 0) and equilibrium (when O = 0) probability densities at T = 7 t/ 2 and is 
expressed as a = ( Vt/h ) sin 2 0. The numerical results for a as a function of 0 are obtained 
by solving Eqs. (5.11) for T\ = 5.56 x 10 _7 s, n = 8, h z = 0 and different values of h and 
Q. The solutions of these equations, i.e., the pairs of angles 9 n = 0{t n ) and ip n = ip(t n ), 
are determined at the moments of time t n = 10 5 + nAt (this choice of the initial time 
guaranties that the transient processes are completed) with n = 1, N, N = 10 11 and At = 
10~ 2 . Finally, the numerical values of the probability densities P s t|^=71-/2 and Pok=7r/2 are 
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FIG. 4. (Color online) Dependence of a on 0 for f l/h = 0.5 (1) and £l/h = 1.5 (2). The solid curves 
correspond to the theoretical result a = (fl//i)sin 2 0 and the symbols represent the numerical 
values of a. The latter are obtained by solving Eqs. (5.11) for h = 2.5 x 10 —2 , fl = 1.25 x 10 -2 (1) 
and = 3.75 x 10~ 2 (2); the other parameters are given in the text. 

calculated as N m \ci^o/(A@A i ^N) and AT m |n = o/(A0A\I/iV), respectively. Here, N m is the 
number of pairs (among total N pairs) satisfying the conditions 9 n G [mA0 — A0,mA0] 
and ip n £ [7r /2 — AT/2,7r/2 + AT/2), in which the parameters m, A0 and Ad* are chosen 
to be m — 1,155 and A0 = AT = 7r/155 (note also that ^ m = N). 

As is illustrated in Fig. 4, our simulation results for the ©-dependence of a are in a very 
good agreement with the theoretical prediction. This leads to the following conclusions. 
First, the solution (5.10) of the Fokker-Planck equation (5.1) correctly describes the long¬ 
time behavior of the rotational motion of nanoparticles in a viscous fluid. Second, since the 
difference (P st — P 0 )|^ =7r / 2 °f the second order in nh, the effective Langevin equations (5.11) 
can be used to predict and study subtle rotational effects. And third, the representation 
Pst = Pst (0, T), which is the key assumption in our analysis, holds not only at fl/h < 1 (as 
it could be expected from the noiseless case), but also at f l/h > 1. The last means that the 
periodic regime of rotation does not influence the steady-state probability density P st . 

The numerical solution of Eqs. (5.11) is then used to determine the average values of 
the nanoparticle magnetization and nanoparticle angular frequency, (/i) and (y). They are 
calculated as (/i) = (1/7) ^- =1 cos 6 1 ; and (x) = 1 - (1/0 Yli=\ where 6 1 * = 9i(t s im ) 

and = 'Oi(tsim) are the polar and lag angles in the i-th run, t sim is the simulation time, 
and l is the total number of runs. In our simulations we set f S i m = 10 6 and l = 10 5 ; the other 
parameters are the same as in Fig. 2. It should be noticed that since t sim is large enough, 
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FIG. 5. (Color online) Average value of the reduced angular frequency of nanoparticles, (\), as a 
function of for h z = 0, h = 1 and different values of the parameter k. The solid curves show the 
results obtained by numerical solution of Eqs. (5.11), and the dashed curve represents Eq. (4.7), 
which corresponds to the noiseless case ((x)U=oo = x)- 

the statistical properties of angles 6 i and ipi do not depend on their initial values. 

Using this approach, we observed that (/r) = 0 for all finite values of the inverse temper¬ 
ature parameter ft. At first sight, this result is in disagreement with the behavior of /i in 
the noiseless case (when ft = oo), see Fig. 2, because (/i) at large ft should approach /i. If 
Ij , is numerically determined for small but non-zero values of h z (e.g., h z = 10 -3 in Fig. 2), 
then (/i) at /c 1 /h z indeed approaches /i. However, since /j is mathematically defined 
as h z — > 0, such an approach is impossible for any finite ft. In fact, the periodic regime of 
nanoparticle rotation, which exists in the noiseless limit at h z = 0 and Q/h > 1, is degener¬ 
ate: lini/, z ^o± The thermal torque of arbitrary strength completely destroys this 

regime, leading to (/i) = 0. 

The influence of thermal torque on the average angular frequency of rotation of nanopar¬ 
ticles driven by a circularly polarized magnetic held (when h z = 0) is illustrated in Fig. 5. 
As seen, the average angular frequency (x) is strongly affected by thermal torque (the less 
the parameter ft, the more the torque strength) and exhibits a remarkable dependence on 
the driving held frequency Q. Since = 0, see Eqs. (5.11), the dependence of (x) 

on ft and is a purely nonlinear effect. Its most striking manifestation is that thermal 
fluctuations can both increase and decrease the angular frequency of nanoparticles as com¬ 
pared with the deterministic case. Specifically, if is the solution of the equation e = 0 
(e = (x)U=oo — (a)) with respect to H, then thermal fluctuations decrease the frequency of 
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FIG. 6. (Color online) Dependence of e on D for different values of k. The simulation parameters 
are chosen as in Fig. 5. 

rotation (e > 0) when fl < f2 c , and increase it (e < 0) when O > (see Fig. 6). Note that 
ff c grows and |min c| decreases as k becomes smaller, and e approaches zero at large fl. 

By solving the effective Langevin equations (5.11) numerically, we investigate the role 
of thermal torque in the nanoparticle dynamics induced by the precessing magnetic held 
(when h z ^ 0). Before we proceed with the analysis of thermal effects, we recall that in 
the noiseless case the steady-state dynamics of nanoparticles has a precessional character 
described by constant polar and lag angles (4.3). As a consequence, in this case y — be., 
the angular frequency of precessional rotation of nanoparticles coincides with the magnetic 
held frequency, and the ^-component of the reduced nanoparticle magnetization is given by 
Eq. (4.4). 

Because Eqs. (5.11) are nonlinear, the thermal torque essentially influences the average 
characteristics of the precessional motion of nanoparticles. In particular, due to its action, 
the average angular frequency of precession becomes less than the magnetic held frequency, 
i.e., (x) < 1 for all finite values of k. Moreover, the numerical simulations show that 
(x) is a monotonically decreasing function of \/k with (x)| 1 / re=0 = 1 and (x)|i/«=oo = 0. 
The average frequency (x) also decreases monotonically with increasing (see Fig. 7), and 
(x) —* 0 as —* oo for each finite n. An important feature of these dependencies is that they 
decrease more slowly with increasing h z . This fact suggests the existence of the characteristic 
frequency Ho, which separates two qualitatively different behaviors of (x) as a function of 
h z . Namely, if ff < then (y) monotonically decreases as h z increases, and (y) exhibits a 
non-monotonic dependence on h z if 11 > fl 0 , as shown in Fig. 8. It is important to emphasise 
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FIG. 7. (Color online) Average value of the reduced angular frequency of precession of nanoparticles 
as a function of 17 for h = 1, k = 5 and different values of the parameter h z . 



FIG. 8. (Color online) Average value of the reduced angular frequency of precession of nanoparticles 
as a function of h z for h = 1, k = 5 and two values of 17. In this case, the characteristic frequency 
17o that separates the regimes of monotonic (17 < 17o) and non-monotonic (17 > 17o) dependence of 
(x) on h z is equal to 0.78. 

that all these remarkable properties of the average frequency of precession of nanoparticles 
result from thermal fluctuations; in the noiseless case (x) = 1. 

Finally, the dependence of the average reduced magnetization (m z ) on 12 and k is il¬ 
lustrated in Fig. 9. As seen, (m z ) approaches the theoretical result (4.4) as n grows, 
and (m z ) almost does not depend on 14 at relatively small k. Since the limit 17 —> oo 
corresponds to the absence of the rotating magnetic held, from Eq. (3.10) one obtains 
(m z ) lo-foo = 27 t cos 0P o (0)d0 = L(kIi z ), where L(x) = cothx — 1/x is the Langevin 
function. In particular, for curves 1, 2 and 3 the function L(nh z ) approximately equals 
0.778, 0.438 and 0.099, respectively. 
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FIG. 9. (Color online) Frequency dependence of the average reduced magnetization for h z = 0.3, 
h = 1 and different values of the parameter k. The dashed curve represents the theoretical result 
(4.4) for the noiseless case. 

VI. CONCLUSIONS 

We have studied both analytically and numerically the rotational properties of ferromag¬ 
netic nanoparticle in a viscous fluid driven by a precessing magnetic held. Our approach 
is based on the system of multiplicative Langevin equations for the polar and azimuthal 
angles of the nanoparticle magnetization frozen into the massless nanoparticle. From these 
equations, approximating the Cartesian components of the random torque by Gaussian white 
noises and interpreting them in an arbitrary way, we have derived the corresponding Fokker- 
Planck equation. By associating the stationary solution of this equation with the Boltzmann 
probability density, we have established that the basic system of Langevin equations should 
be interpreted in the Stratonovich sense. Within this framework, we have reproduced the 
known system of effective Langevin equations, which is simpler than the basic one, and have 
shown that the statistical properties of its solution do not depend on the interpretation of 
multiplicative white noises. 

Using the system of effective Langevin equations and the corresponding Fokker-Planck 
equation, we have calculated the average angular frequency of precession of nanoparticles 
and the average magnetization of nanoparticles in the z-direction, and have analyzed their 
dependence on the model parameters. In the noiseless limit, the dependence of these quanti¬ 
ties on the rotating held frequency and amplitude is different whether a constant component 
of the magnetic held is zero or not. In the former case, it has been shown both analytically 
and numerically that the angular frequency and magnetization depend only on the ratio 
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of the rotating field frequency to the rotating held amplitude and, starting from a certain 
value of this ratio, these dependencies become strongly nonlinear. The most remarkable 
property of the above mentioned quantities is that their sum is strictly equal to 1 (in di¬ 
mensionless units) for any rotating held. In the latter case, when the steady-state rotation 
of nanoparticles has the precessional character, we have derived a general expression for the 
nanoparticle magnetization and have observed that the frequency of nanoparticle precession 
always coincides with the rotating held frequency. 

The influence of thermal fluctuations on the rotational dynamics of nanoparticles is in¬ 
vestigated by numerical integration of the system of effective Langevin equations. To verify 
these equations, we hrst calculated the difference between the steady-state and equilibrium 
probability densities of the nanoparticle orientation, which arises from a slowly rotating 
magnetic held of small amplitude. Then, by comparing the numerical results for this differ¬ 
ence with the results obtained from the analytical solution of the Fokker-Planck equation, we 
have confirmed the validity of effective Langevin equations. Finally, using these equations, 
we have observed a number of interesting thermal effects. In particular, it has turned out 
that the deterministic regime of rotation of nanoparticles, which exists when a constant mag¬ 
netic held is infinitesimally small and the rotating held frequency exceeds the critical one, is 
completely destroyed by thermal fluctuations. But the most important observation is that 
thermal fluctuations can play a constructive role in the precessional dynamics of nanoparti¬ 
cles. The non-monotonic behavior of the average angular frequency of nanoparticle rotation 
as a function of the constant magnetic held strength supports this statement. 
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